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Abstract 

We calculate the vacuum energy of a spinor field in the background 
of a Nielsen-Olesen vortex. We use the method of representing the vac- 
uum energy in terms of the Jost function on the imaginary momentum 
axis. Renormalization is carried out using the heat kernel expansion and 
zeta functional regularization. With this method well convergent sums and 
integrals emerge which allow for an efficient numerical calculation of the 
vacuum energy in the given case where the background is not known ana- 
lytically but only numerically. The vacuum energy is calculated for several 
choices of the parameters and it turns out to give small corrections to the 
classical energy. 



1 Introduction 

Quantum corrections to classical background configurations are a topic of con- 
tinuing interest. At present it is stimulated by the observation made in lattice 
calculations that the field configurations responsible for confinement are domi- 
nated by monopole or string like configurations. Another motivation comes from 
the stability analysis of Z-strings with respect to fermionic fluctuations. 

During the last years quantum corrections to string like configurations have 
been investigated quite actively. In pQ for the background of a finite radius mag- 
netic flux tube in QED the vacuum energy of a spinor was calculated. There 
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the method of representing the vacuum energy in terms of the Jost function of 
the related scattering problem taken on the imaginary momentum axis was ap- 
plied. This method has been developed earlier in [2j for spherically symmetric 
scalar background fields. The specific example considered in pQ was a homoge- 
nous magnetic field inside the flux tube. This investigation was extended to more 
complicated profiles of the magnetic field inside the flux tube in [3J Hj. In [S] a 
magnetic background was considered which depends only on one spatial coor- 
dinate. In addition, this dependence is of a form that the corresponding wave 
equation has an explicit solution so that the vacuum energy could be calculated 
quite easily. A more general approach to fermionic vacuum energy was taken in 
[S| where general bounds on the fermionic determinant were obtained. Another 
approach was taken in [7j where several profiles of the magnetic background were 
considered and compared with the derivative expansion. In jH] the limit of a 
strong magnetic field was investigated in more detail. 

There is still an interest in backgrounds of infinitely thin strings which consti- 
tute a singular background. Typical for these configurations is the need to apply 
the method of self adjoint extensions. The object to consider in such examples 
is the vacuum energy density per unit volume rather than the 'global' vacuum 
energy which is for a string the density per unit length. For recent investigations 
see, for example, [Oj and The problem with infinitely thin strings is that 

their classical energy is infinite. Also, there are additional counter terms. In the 
language of heat kernel expansion there are additional contributions to the heat 
kernel coefficients, for instance coefficients with half integer number, which reside 
on the surface where the singularity is located. Another related example was con- 
sidered in |3j where the background is a finite radius flux tube with the magnetic 
field concentrated on the surface of the string. The advantage of singular (and 
non smooth like in pQ) backgrounds can be seen in the usually quite explicite for- 
mulas for the quantum fluctuations. So in EH E] only Bessel functions appear 
and in m E] hyper geometric functions in addition. 

In general, physical background configurations should have a finite energy, 
hence strings should have a non zero radius. A typical example is the Nielsen- 
Olesen string [TI]. But here, not only the vacuum fluctuations have to be calcu- 
lated numerically, but even the background itself. The problem appears to have a 
calculational scheme which does not need explicite formulas and which allows for 
efficient numerical evaluation. There the main problem comes from the ultravio- 
let divergencies. In analytical terms it is well known how to handle them. First 
one has to introduce some intermediate regularization. After that one has to sub- 
tract the counter terms and finally to remove the regularization resulting now in 
a finite result. However, consider the last step in zeta functional regularization. 
Here one has to perform an analytical continuation. Or let us consider some 
cut-off regularization, where one has to remove the contributions proportional to 
non negative powers of the cut-off parameter. In principle such a procedure can 
be done numerically (there are some examples) but this is quite complicated and 
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ineffective. It is better to transform the expression for the vacuum energy in a 
way that the final removal of the regularization can be performed analytically 
so that only well convergent sums and integrals remain. Such a method had 
been developed in |T21 E] and in [TJ H3j applied to strings of finite radius. The 
method is based on a representation of the regularized vacuum energy in terms 
of the Jost function of the related scattering problem taken on the imaginary 
momentum axis. Another method, using phase shifts and momenta on the real 
axis was used in (see also ^3] and references therein), mainly for spherically 
symmetric backgrounds. Also, using phase shifts, the case of a color magnetic 
vortex was considered in In a similar way in the vacuum energy for an 
electroweak string had been considered, where, however, a step profile was taken 
in the final stage. 

A completely different method is worth to be mentioned. In ^HJ E|] world 
line methods were applied to the calculation of vacuum energy which have the 
advantage not to rely on separation of variables and therefore to be applicable 
for much more general background configurations. 

In the present paper we calculate the vacuum energy of a fermion in the 
background of a Nielsen-Olesen string. We use the method of representing the 
regularized vacuum energy in terms of the Jost function on the imaginary mo- 
mentum axis which was developed in For the Nielsen-Olesen vortex the 
background potential is given only as a numerical solution of the corresponding 
equations of motion. We use zeta functional regularization and determine the 
counter terms from standard heat kernel expansion. The renormalized vacuum 
energy is divided into two parts, the 'finite' and the asymptotic ones, by sub- 
traction and addition of some first terms of the uniform asymptotic expansion of 
the Jost function which is obtained using the Lipmann-Schwinger equation. In 
the asymptotic part the analytic continuation in the regularization parameter is 
performed analytically and well convergent double integrals remain. The 'finite' 
part is represented as a well convergent sum and integral involving the Jost func- 
tion which is calculated from the numerical solutions of the corresponding wave 
equation. Using these tools, the dependence on the parameters of the considered 
model is investigated numerically. 

The paper is organized as follows. In the next section the basic notations 
for the considered model are introduced, the renormalization is discussed and the 
basic formulas for the representation of the vacuum energy are given. In the third 
section the asymptotic part of the Jost function is derived and the asymptotic 
part of the vacuum energy as well. In the fourth section the 'finite' part of the 
vacuum energy is derived and the convergence properties are discussed. Sect. 5 
contains the numerical part of the work and Sect. 6 the conclusions. 
We use units where h = c = 1. 
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2 Basic notations 



The Abelian Higgs model contains a f/(l) gauge field, and a complex 

scalar field, The action is 

^ = /rf 4 x I - A (W - f )1 , (1) 

where _D M = <9 M — igA^ is the covariant derivative and F Ml , = d^A v — d u A^ is the 
field strength, for more details see [20] • The vacuum configuration is given by 
= 0, $ = ?ye lc / a/2 where c is some constant and r\ is the Higgs condensate. 
This configuration has zero energy. A configuration with finite non zero energy 
must be at spatial infinity in the vacuum manifold. Hence asymptotically the 
gauge potential must be a pure gauge and the scalar field must tend to a constant 
times an angular dependent phase. A configuration of such type is the Nielsen- 
Olesen string [TT]. In cylindrical coordinates (x,y,z) — > (r,(p,z) one makes the 
ansatz 

$ = -^/(r)e^, qA p = nv{r), (2) 

where A^ is the angular component of the vector potential. The profile functions 
satisfy the boundary conditions 

/(0) = v (0) = 0, f(r) -> l,v(r) -> 1 as r -> oo. (3) 

Here n is the winding number and gives at once the magnetic flux in the string. 
In the following we consider n = 1 only. The equations of motion imply 

f^r) + ~f(r)-±(l-v(r)) 2 f(r)+\r) 2 (l-f(r) 2 )f(r) = 0, 

V "( T )- \ v >{r) +q Yf{r) 2 {l-v{r)) = 0. (4) 
r 

After a rescaling, r = p/(qrj), these equations depend in fact only on the com- 
bination (3 = 2X/q 2 , which is at once the squared ratio of the Higgs and vector 
masses. For (3 = 1 the system exhibits an additional symmetry and can be re- 
duced to two first order equations (Bogomolny equations). For (3 < 1 the system 
is stable for all values of the flux n = 1,2,..., for j3 > 1 only n = 1 is stable (a 
vortex with n > 1 decays into n vortexes with n = 1). 

The system of equations (jlj) together with the boundary conditions Q does 
not have an analytical solution and one is left with numerical methods. The 
simplest way is to set the derivatives in zero, i.e. the constants a and b the 
expansion for small r, f(r) = ar + . . . and v(r) = br 2 + . . ., and to numerically 
integrate the equations from r = to larger values of r. One has to adjust these 
constants in a way that for large r the asymptotic values / = 1 and v — 1 are 
approached. Examples are shown in Fig. 1 for two values of (3 with qrj = 1. The 
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Figure 1: The profile functions of the Nielsen- Olesen vortex as functions of the 
radius r for qr] — 1. 

corresponding values of the derivatives are (a, b) = (0.26817, 0.17481) for ft = 0.3 
and (a, b) = (0.79958, 0.42848) for (3 = 6. As the asymptotic values for r — > oo 
are reached with exponential speed it is in fact sufficient to consider r up to ~ 15. 
sufficient to consider r up to ~ 15. 

The classical energy (more exactly, the energy density per unit length) of 
these configurations is given by 

^ciass = n Tdrr (^'(p) 2 + V 2 (f{pf + 1(1 - v(p)ff(pf) 
Jo \q z r z \ r z J 

+^ 4 (1 - f(p) 2 ) 2 ) • (5) 

We calculate the background for qrj = 1 and restore later the general setting by 
inserting p = qrjr in (JHJ). 

The spinor is taken as a four component Dirac spinor with a coupling to the 
background given by the Lagrange density 

L = -{^(t^D,-f e \<t>\)^. (6) 

This model is chosen for the reason of simplicity. It provides a coupling of the 
spinor not only to the vector but also to the scalar background and it is motivated 
by the Yukawa coupling in the fermionic sector of the standard model. We note 
that the interaction in © is gauge invariant and that the coupling constant f e 
is dimensionless. However, this model has a drawback. Since there is only one 
spinor we are forced to take the absolute value, | $ |= J (5R$) 2 + ($5$) 2 , of the 
complex scalar field (in the standard model there we don't need to do so). As 
a consequence, the interaction is non polynomial. This is a complication, for 
example if considering the vacuum energy together with the dynamics of the 
background. Also, there is a need for an additional counter term (see below). 
However, because we are not going to consider the dynamics of the background, 
this complication does not affect the calculation of the vacuum energy if we let 
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enter the additional counter terms into the classical energy with a coefficient 
which is put equal to zero after the renormalization is carried out. In addition, if 
we consider the same problem of calculating the vacuum energy in the standard 
model then there is no such complication due to its renormalizability. 

The background is static and due to the translational invariance in direction 
along the z-axis and the rotational invariance around the z-axis (m = 0, ±1, ... is 
the orbital angular momentum), the corresponding momenta can be separated, 
^ — > e ~ l Po^°+W3x 3 f^ie ^^Jj^j, After that the Dirac equation decouples and one 
of the resulting pair of two component equations may be represented as 



po - M r ) 



a 



V -3F r > P0 + K r ) J V W2 / 

The other pair is obtained by changing the sign of \x. Here the notation /i(r) = 
f e -j=f{r) has been introduced. With respect to the spinor this is a radius de- 
pendent mass density. From its value at infinity we define the spinor mass, 
m e = / e ^|- For a constant mass density the problem is the same as for a pure 
magnetic flux tube which appears as a special case in this way. 

The vacuum energy of the spinor is given by the general formula 

Equant = ~ ~ E e \n) S - (8) 
Z (n) 

Here s is the zeta functional regularization parameter and e< n \ are the one particle 
energies. For simplicity we dropped the arbitrary constant /i which is usually in- 
troduced to adjust the dimension in zeta functional regularization. The quantum 
numbers (n) include the sign of the one particle energies (all enter with positive 
sign), the spin s z (two projections which we can account for by summing over the 
sign of /i(r)), the orbital angular quantum number m and the radial quantum 
number n r (assuming for a moment the system being inserted into a large cylin- 
der). In addition there is the momentum ^3 which can be integrated according 
to 

" ?P3 ( 2 2\4— = 2(1-) r ( g ~ 1) = x 2{X- S )£s_ 

00 2tt v 3 J 2^T(s - \) Aits 

with C s = 1 + s(2 In 2 — 1) + . . . and we arrive at 

Equant = E (e SMr ) 2il ~ S) ■ (9) 

In order to get rid of the large cylinder we proceed as in 0. We consider the cylin- 
drical scattering problem associated with Eq. (J7j) and define the Jost functions, 
fm{k), with po = sjvn? e + k 2 . We rewrite the sum over n r by an integral, tend the 
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large cylinder to infinity and after dropping the Minkowski space contribution we 
end up with 



^uant = /E E / dk (k 2 - mfj — In f m {ik) (10) 



which is our expression for the regularized ground state energy. Here the inte- 
gration is turned to the imaginary axis, more specifically, it goes along the cut 
resulting from (k 2 — m 2 e ) 1 s . In order to fully exploit the symmetry we renumber 
the orbital momenta by v according to 

/ ""5 (m - 0) nu 
m = < i . , (11) 

\ -v - \ (m < 0) y J 

with v = |, |, . . . and two signs of the orbital momentum are to be taken into 
account. After rotation k — > ik we define p = \Jk 2 — m 2 and for positive orbital 
momentum Eq. (J7J) can be rewritten in the form 



ip — fj,(r) 



d_ _ u-l/2-v(r) 



y ! %+Mr) JUJ"°- (12) 

It can be seen (see below, Eqs. ()44j) and (|45jl) that for the Jost function a change 
in the sign of \x corresponds to complex conjugation and a change in the sign of 
the magnetic background, v(r) — ► — f(r), corresponds to an exchange of positive 
and negative orbital momenta. As a consequence, in Eq. (fTU|) we have a factor 
of 4 and we take the real part of the half sum of positive and negative orbital 
momenta. 

The next step is the renormalization. We use the standard heat kernel expan- 
sion according to which the divergent part of the ground state energy is given by 
(see, e.g., [T] and we drop a ) 

^div = [ I + in _L - 1 ] at l — ( - + In 4r - 2 ) a 2 , (13) 

32tt 2 \s ml J 32tt 2 \s m 2 e J 2 ' v ; 

where a« are the standard heat kernel coefficients. We define the renormalized 
ground state energy by s — > 

Ex e xi Equant -£/ (^-4) 

in the limit s — > 0. The definition (fT4"|) of the renormalized ground state energy 
is chosen in a way that E ven vanishes if the spinor mass m e taken as a parameter 
independent from the background becomes large. This is the so called 'large 
mass' normalization condition. It has been discussed in detail in [221 E3| and in 
[2*4*] it was shown to be equivalent to the known 'no tadpole' condition. 
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The heat kernel coefficients can be calculated using known methods from the 
squared Dirac operator, 

Mr)) (ip- A*(r)) = —D 2 + - i/.V^; I $ I ~fl I $ | 2 

= -D 2 + K (15) 

with = |[7 M , 7"]. The general expressions for the relevant heat kernel coeffi- 
cients are 



ai = Tr J d 2 x (-V), 



«2 = Tr I d 2 * (-^ + V - Iav) , (16) 

where the trace is over the gamma matrices. Inserting for V the result is 

00 

ai = — 87r ^ dr r (ji{r) 2 — m 2 j , 

^^ + ^(^V) 2 + (Mr) 2 -m 2 ) 2 )) (17) 

(note that we consider densities per unit length of the string). From Eq. (J 13)) we 
obtain the divergent part of the vacuum energy in the form 

£div = 1 - 2 g (ln(m e /2) + 1) J dr r ( IvVf + Sjtf {r f 




Airs J \ 3 r 



~^J drr 4 (/W 2 "l)+0( S ). (18) 

n 



The interpretation is as follows. From the f'(r) 2 -term we have the standard 
renormalization of the electric charge q in the classical action (|5f . A renormaliza- 
tion of the scalar coupling A absorbs the divergence proportional to fir) 4 . After 
that the remaining freedom is in a change of the condensate rj. However, this is 
obviously insufficient to absorb the remaining parts of S dlv . In this way, in the 
classical energy an additional structure must be present. As suggested from Eq. 

(fi"5|) . it is necessary to introduce at last a term proportional to \ $ |) into 
the action (JTJ) and into the classical energy (JEJ) as well. It should be noted that 
such a term is gauge invariant and that it has the correct dimension. However, 
it represents another non polynomial interaction. This is not surprising as the 
model given by Eq. (jUJ) itself contains a non polynomial interaction. Accepting 
this we can finish the renormalization if we put the coefficient in front of this 
term equal to zero after performing the renormalization. 
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In order to perform the limit s — > in Eq. (|T4*|) we need to rewrite the 
regularized ground state energy. For this we define the asymptotic Jost function 
as the part of its uniform asymptotic expansion for v — * oo, z = - fixed, which 
includes all powers up to z/~ 3 , 

In /,,(<*)= In + (19) 
Using In fjf we divide the energy into 

E mn = S { + (20) 

with the 'finite' part, 

= - E / dk (k 2 -ml — (In -In /"(<*)) (21) 

and the 'asymptotic' part, 

£T=i E Hdk (tf-ml)^ ±lnfr(ik) - E«*. (22) 

7T Jm e v 'OK 

u ~ 2 >2 

The sum and the integral in £ f , Eq. (|21|). are finite by construction of the asymp- 
totic Jost function. Therefore we could put s = therein. In the asymptotic part, 
i? 88 , Eq. (|22|) . the analytic continuation in s has to be still performed which will 
be done in the next section. 



3 The asymptotic part of the vacuum energy 

The asymptotic part of the Jost function can be derived from the Lipmann- 
Schwinger equation just generalizing the procedure developed in pQ. The operator 
AV there in Eq. (28) reads now 



/i(r) - — 



AP= ( ' ]. (23) 



Again, we have to perform iterations up to the fourth order in the operator AV. 
Using the formulas given in £Q one arrives at the representation 



y j 3 3n p 

wrw = /yEE^A (24) 

n=l j=n 
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with t = —, — = . The coefficents are 

X n = I w ( r ) 2 + ir 2 (^(r) 2 -m e 2 ), 

■\r 1 / \ 2 1 9 / / \ 2 1 / \4 1 o / \ 2 / 2 2\ 

A33 = 7u(r) - -r f (r) - -v{r) - -r v(r) [ji - m e ) 

4 o o 4 

-^ 4 (/i'(r) 2 + (/i 2 -m 2 ) 2 ), 

a/- / \ 2 "^2//\2 3 / \ 4 3 n / \Q/ 2 2\ ^2/2 2\ 

^37 = f^) 2 -^(r) 4 + ^r 2 (^ 2 -m 2 ), 

^39 = -y|^(^) 2 - (25) 

The dependence on v(r) is the same as in 1] (where it was denoted by a(r)), the 
dependence on /x(r) is new. In Eq. (I24|) it had been integrated by parts in order 
to get the shortest representation for the coefficients X n j. 

In £ as , Eq. (|22JI. the analytic continuation in s can be performed by rewriting 
the sum over the orbital momenta v by integrals using the Abel-Plana formula in 
the form as given in the Appendix, Eq. (|53|). From the first part, i.e., from the 
direct integral over u, we get a contribution which just cancels E dlv in Eq. (|22|) . 
So we are left with the contribution from the second part. Here we integrate over 
k using the simple formula and obtain 



with 



CO 



' r' ./ I • e 2 ™ i 
o 



. (27) 



The difference in the right hand side results from the deformation of the integra- 
tion contour in the Abel-Plana formula which gets tight around the cut. There- 
fore the integration starts effectively from v — x. This formula provides the best 
representation to perform the analytic continuation in s. For (n,j) = (1,1) we 
simply note 



E u (x) = — / du - — vV-x 2 . (28) 

v ' x 2 Jx 1 + e 2 ™ v ' 

Also for (n,j) = (1,3) we can put s = immediately. However, we prefer to 
integrate by parts in order to get a representation which is in line with Eq. 
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In a similar way, integrating by parts, we can proceed in all other contributions. 
So the final form is 



2 roo 

E n j(x) = -o / dv fnj{v) \Jv 2 -X 1 



(29) 



with 



/33(f) 
/37(f) 



1 

1+e 2 ™ ' 

( 1 



„(i+ e 2™) y , 



3 fi> fi> (l+e 2 ™) J ) 



/13(f) 
/35(f) 

/39(f) 



1+e 2 ™ y ' 



1 

v v 1+e 2 ™ y / ' 



1 f 1 / 1 / 1 



15 I v \ v \ v \ 1+e 



Using these formulas in Eq. ()26j) we can put s = there and insert for the 
coefficients X n j using Eqs. (j25)l . After rearranging contributions we arrive at 



S 3 



'2 r<x> (if / 

— / — (v [r) 2 hi{rm e ) + r 2 v'(r) 2 h2(rm e ) + v (r) 4 /i 3 (rm e ) 
w(r) 2 (/i(r) 2 — m 2 )/i 4 (rm e ) + r 2 (/i(r) 2 — m 2 )/i 5 (rm e ) 



with 



+r 4 (/i'(r) 2 + (/x(r) 2 - m 2 ) 2 )/i 6 (rm e ) 
= / dv /j(f) V v 2 — x 



(30) 
(31) 



and 



/i(f 
/ 2 (f 
/ 3 (f 
/ 4 (f 
/s(f 
/e(f 



-/11(f) - /13(f) + ^ /33(f) - y /35(f) + ^ /37(f) - ^/sg(f), 

-^/ss(f) + ^/ss(f), 

--/33(f) + ^/ss(f) - jMf), 
1 1 

-g/aaW + ^(f), 

-/11(f) - ^/ss(f) + ^ /37(f), 
-^/ss(f)- 



(32) 



According to these functions we divide the asymptotic part of the ground state 
energy, 

6 

£a * = H £ F> ( 33 ) 
3=1 

into parts which have the meaning, for instance, of the asymptotic part of the 
vacuum energy resulting from the magnetic background, £f s , £f s and S^ s , or 
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from the mixed contribution, £|" s , etc. We finish this section with the remark, 
that the integrals in Eq. (|3U|) are all finite for the considered background of 
a Nielsen- Olesen vortex due to corresponding properties of the functions hi(x), 
their behavior at x — > for instance. 



4 The 'finite' part of the vacuum energy and 
numerical results 

In order to calculate the finite part of the vacum energy defined by Eq. (|2"Tf 
we have, first of all, to set up a numerical scheme for the calculation of the Jost 
function. The Jost function is defined as the coefficients in the asymptotics for 
large radius in the so called regular solution of the wave equation (JZJ), 

* {±) (r) ~ ~ {f^ik^ikr) +7? ) (*)*gi)(*r)} , (34) 



where 



(35) 



(po — y m e + are the two li near independent solutions at r — > oo. The signs 
(±) correspond to the sign of the orbital momentum according to Eq. (jllj) and 
H^' 2 \z) are Hankel functions. In general, 5 is the dimensionless value of the 
magnetic flux, i.e. it is v (r — > oo), which is equal to one in our case. 

From Eq. (|34p. the Jost function can be expressed in terms of the solutions, 



(36) 

where ipi^{kr) are the upper and lower components of \l/(r). Strictly speaking, 
as given by Eq. (J3HJ), the functions fj^(k) depend on r and only for r — > oo 
they tend to the Jost functions. However, for sufficiently large r (larger than 
the scale of the background, say r > 15 in the examples shown in Fig. 1), 
Eq. ({SEt provides a good approximation. In the considered case the background 
approaches its asymptotic values with exponential speed. In the same way the 
difference between the Jost function and the approximation (|3*6^1 is small. 

The regular solutions \&(r) of Eq. (j7J) are defined as becoming proportional 
to the free solutions for r — ► and can be expressed in terms of Bessel function, 

y/po + m J„±i(qr) 
- - , / j 2 ( J <•><> 
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with q = \Jk 2 — ml + m\. We introduced for a moment mo = /i(0) which cor- 
responds to a more general case. In all examples considered below we will have 
m = 0. Note also the factor in front which provides the correct normalization. 

In order to actually carry out a numerical integration of Eq. (J7J) it is useful to 
change to functions with regular initial values at r = 0. At once, we must change 
to imaginary momenta. So we substitute k —>■ ik and we make the ansatz 

, (+) = r| Wp=*sM +) \ ¥ - ) = l J Jv=SSi4P \ m 

+ im 02 +) / V ~WP + im o 02 +) / 



The equations for are 



Or 



-(P + V(r))e« , £ + ^±^ \ w _ 
f-^^ , -(p-WJe-ar _0 (40) 



and 



with e = y(p + W7i e )/ (p — zm ). Further we substitute 
and 





and they obey the equations 

-(P + V*(0)5gre* > f + ^V 1 

l_ 2^+1 

<9r 



4 + , -ip-i^r)) 



-(p-z>(r)) 




(39) 



(41) 



(42) 



The boundary conditions at r = for these functions are 

9{fl)(r)\r= = 1 (43) 



(44) 



i-¥ \ -(p-w)a--" 11 " 1_0 ' (45) 

The above mentioned symmetries can be seen in this representation explicitely. 
The Jost functions read now 

i 



fl + \ik) = ^fij (^i^'W^i-iW + w *9i + \r)K v+ i_ s (kr) 



(46) 
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and 



2 (47) 
with w = \j {p + im e ) (p — imo) and modified Bessel functions. Eqs. f)44j) and 
(145)) together with the boundary conditions (|4"3|) can be solved easily numerically. 
We used the package NDSolve in Mathematica. 

Next we need the asymptotic part of the Jost function. We use Eq. with 
the coefficients X n j, Eq. (|25j). The integrals over r are convergent and can be 
calculated numerically without any problems. Taking into account the mentioned 
symmetries we drop all contributions which are odd under v — > —v and we define 
the subtracted logarithm of the Jost function as 

ln/ sub = ^ (ln/ (+) {tk ) +l nf (-)(ik)) - In/- (48) 

This expression we insert into £ f , Eq. (|21|). After integration by parts we arrive 
at 

£'= E Z (49) 
with the contributions from the individual orbital momenta, 

roo 

Sl= dk El(k), (50) 

and the integrands 

Sl(k) = k ln/ sub . (51) 

7T 

As illustration of the convergence we show £ l u {k) enhanced by a factor of k 3 u 3 in 
Fig. 2 for several first v. It is seen that all curves approach the same limiting value 
which is essentially given by the next term in the uniform asymptotic expansion 
of the logarithm of the Jost function after that included into ln/ as , Eq. (fTT)|) 
taking into account that the order \ is absent for symmetry reasons and the 
order -g is the next one. In this way the integrals over k are fast convergent. 

In Fig. 3 we show the contributions from the individual orbital momenta. 
The corresponding sum in Eq. (|4T?|) is also fast convergent and it is sufficient to 
take a quite small number of the first orbital momenta. 

We would like to add a note on bound states. They appear for imaginary 
momenta, k = in, in Eq. (J7J), in the region where po = \jm 2 e — n 2 is real, 
< po < m e . For such momenta the representation (j3T?j) for the Jost function 
can be rewritten in the form 

= -kr C^G-*) {V-Po + m e y 1 (r)K l/T( i +S) (Kr) 

+ ^/p Q + m e m 2 (r)K v±c _ +s) (Kr)} . (52) 



14 




Figure 2: The integrand £l(k) in Eq. (J50|) multiplied by fc 3 z/ 3 as a function of k 
for several first values of v. In this figure and in the next two the parameters are 
(3 = 0.3, q = 0.5, f e = 1, rj = 1. 



The functions ^2 are solutions of the Eqs. (JJJ) with the initial conditions (J37j) . 
As both, equation and initial conditions are real for the considered momenta the 
solutions are also real. In this way the expression in the figure brackets in Eq. 
(152)1 are real. Their zeros just determine the location of the bound states. As an 
example we plot in Fig. 4 these functions for the two lowest orbital momenta. For 
orbital momentum m = we see one bound state, for negative, m = —1, none. 
The appearance of the bound states can be explained easily. Without scalar 
potential, i.e., for a constant /x(r) = m e the spinor moves in a pure magnetic field 
and in the state where the magnetic moment is antiparallel to the magnetic field 
its coupling just compensates the lowest Landau level resulting in a zero mode 
known from [21]. Switching on the scalar potential the spinor feels an additional 
attraction and becomes a bound state. As concerns the vacuum energy these 
bound states are outside of the integration region in Eq. (jTU|l and therefore 
they do not complicate the calculations. As shown in j2j they are accounted for 
automatically in the representation of the vacuum energy in terms of the Jost 
function. The only we have to take care of is a possible appearance of bound 
states on the imaginary fc-axis in Eq. (JTU)) above m e . This would require a zero 
of the Jost functions in representation ()46j) or ()47)1 for real k there. But these 
are genuine complex expressions and it can be shown that they do not have such 
zeros. 
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Figure 3: The contributions from the individual orbital momenta El to E l in Eq. 

flU. 
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p = 


0.3 


0.47913 


1.50653 


1.91653 


3.90220 


p = 


6 


0.97224 


2.64606 


3.88896 


7.50726 



Table 1: Two examples for the classical energy and its constituent parts for 
q = 0.5, i] = 1. 

5 Numerical results 

We investigated numerically the vacuum energy for values of (3 ranging from 
(3 = 0.3 to (3 = 6 for some choices of the remaining parameters. The convergence 
properties of the sums and integrals involved have been discussed in the previous 
section. Here, let us first discuss the weight of the individual contributions. 

We start with the classical energy which is given by Eq. (JHJ). It can be divided 
into three parts. The first one, -E m ag, is the energy of the magnetic field. The 
second one, E cov , results from the covariant derivative in Eq. ((TJ and is given 
by the second and the third contributions in (jSj). The third one, E sca \, is the self 
energy of the Higgs and is given by the fourth contributions in (jSJ). We show in 
Table 1 two examples, one for smallest coupling of the scalar, A = /3q 2 /2 and the 
other for the largest value considered. In both cases as well as in all intermediate 
ones, the scalar self energy contribution is larger than the magnetic energy and 
than E cov . 

Next we consider for the parameters q = 0.5, f e = 1, r\ = 1 the constituent 
parts of the asymptotic part of the vacuum energy, i.e., the £^ as defined in Eq. 
f!33|) . These quantities are shown in Table 2. As it can be seen the contributions 
£g s and £g s are dominating. They result from the scalar background, see Eq. ([30)1 . 
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Figure 4: The logarithm of the expressions in the figure brackets in Eq. (152)) 
as function of k for lowest orbital momenta. The peak indicates a bound state. 
Its location is k = 0.039232. 
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-0.0017591 


-0.0061022 




ft = 


6 


-0.003898 


-0.002483 


-0.006208 





Table 2: the constituent parts of the asymptotic part of the vacuum energy for 
q = 0.5, / e = 1, 7] = 1. 



But generally, all have at last a numerical smallness of two orders of magnitude. 

Now, in Table 3, we represent the result of the calculations of the complete 
vacuum energy for all considered examples. For two cases we also represent the 
parts of the vacuum energy as a function of ft, see Figs. 5 and 6. From these 
results it can be seen that there seems to be no general rule which part of the 
vacuum energy is dominating and which not. Moreover, in dependence on the 
parameters one part may be larger than other and smaller as well. It is even 
impossible to say something definite about the sign of the vacuum energy, it may 
change although in most cases it is negative. 
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Table 3: The constituent parts of of the vacuum energy for all considered exam- 
ples. 
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Figure 5: The vacuum energy as a function of (3 for the q = 0.5, f e = 1, r] = 1. 
In order to represent all quantities within one plot the classical energy is divided 
by 200. 

6 Conclusions 

In the present paper we calculated numerically the vacuum energy of a fermion 
in the background of a Nielsen-Olesen string. As for renormalization we used 
standard zeta functional regularization and determined the counter terms from 
the first heat kernel coefficients (up to t^). It turned out that from the renor- 
malization a counter term appears which is not present in the initial action. It 
is gauge invariant and it has the correct dimension but it represents a non poly- 
nomial interaction. This has to be considered together with the non polynomial 
interaction present in the model itself as discussed in Sect. 2. 

The numerical investigations have been performed using methods developed 
in the papers 00]. In the present paper the background is given purely numer- 
ically in difference to the previous papers where it had been given analytically, 
for example as a step function or some Gaussian profile. It has been demon- 
strated that the computational scheme used here is well suited to work with such 
backgrounds. This is a step forward to physically really interesting problems. 
In the considered model the stability of the background is given by topological 
arguments and for a realistic choice of the parameters the quantum contribution 
is small. This smallness has two sources. The one is the smallness of the cou- 
pling constants which appears in front of the quantum contribution relative to 
the classical one. The second one is a purely numerical smallness of about two 
orders of magnitude as discussed in Sect. 5. It is present even if the parameters 
and couplings are all of order one. It is obviously connected with the dimension 



19 



ren 



Figure 6: The same as in figure 5 but for q = 2. 

and the ultraviolet renormalization. So in |12j it was found for a one dimensional 
background that in (1+1) dimension the vacuum energy is by one order of magni- 
tude larger than for the same problem in (3+1) dimensions. Comparing formulas 
(29) and (32) in [T2| the interpretation is suggested that this difference is due to 
the one additional ultraviolet subtraction in the (3+1) dimensional case. 

As seen in Sect. 5, for the dependence of the vacuum energy on the various 
parameters no general rule can be seen so far. Even the sign of the vacuum energy 
changes. The same applies to the relative weight of the individual contributions. 
So sometimes the asymptotic part is dominating, in other cases, however, the 
'finite' part is larger. In general, they are of the same order. From this one can 
conclude only that in the given background there is no small parameter which 
could allow for some approximative scheme. So for instance, if the asymptotic 
part is dominating one could hope to get a good approximation by including 
higher orders of the uniform asymptotic expansion of the Jost function into the 
asymptotic part of the vacuum energy and neglect the 'finite' part of it which is 
the numerically much harder part of the problem. The examples in ^3] have been 
of a kind suggesting this way in contrast to the example in the present paper. 

For the considered model of a spinor in the background of the Nielsen-Olesen 
vortex, due to its smallness, the vacuum energy has only a very small influence on 
the dynamics of the background. Hence in the considered case the vacuum energy 
is of limited physical importance. However, its calculation gave new insights into 
the structure of such calculations and demonstrated the power of the methods 
used. A next step could be to apply them to the Z and electroweak strings where 
the stability is not guaranteed by topological arguments and where the stability 
issue is not finally settled with respect to the fermion contributions [2lH I17j. 
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Appendix 

The Abel-Plana formula used in Sec. 2 reads 



£ /M-f^H + / 1+ t,./ W -/'"' J> < 53 > 



--hi- o o 



The following formulas are used in the text, 



oc 



f dH e - - -^ £(^±1^- gz__ (54) 



HIT 



with t = 1/y 1 + (kr/u) 2 . They can be easily derived, see also (C3) and (C2) in 

Hi- 
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